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PROBING THE EXPLOSION MECHANISM OF SUPERNOVAE 
BY RADIOACTIVE DECAY GAMMA- AND X-RAYS 

K. Maeda 1,2 
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ABSTRACT 

Gamma- and X-rays resulting from radioactive decays provide a potentially powerful tool to investi- 
gate the explosion physics of supernovae, since the distribution and the amount of radioactive isotopes 
arc strongly dependent on properties of the explosion. In this paper, expected features of these high 
energy emissions are presented for bipolar jet-induced explosion models, which are recently favored for 
very energetic supernovae and even for typical Type Ib/c supernovae. It is shown that combination of 
various observations, i.e., line-to-continuum ratio, photoelectric absorption cut-off energy, line profiles 
and luminosities, allows the unique determination of the explosion energy, the amount of radioactive 
56 Ni, the explosion geometry, and even the viewing orientation. 

Subject headings: gamma-rays and X-rays - radiative transfer 



1. INTRODUCTION 

At the supernova (SN) explosion, various heavy el- 
ements/isotopes are produced as a strong shock wave 
passes through the inner region of the progenitor star. 
Among these newly synthesized products, important is 
radioactive 56 Ni (Clayton et al. 1969). The decay chain 
56 Ni — > 56 Co — > 56 Fe produces 7-ray lines with aver- 
age energy per decay ~ 1.7 MeV ( 56 Ni decay with an 
e-folding time of 8.8 days) or ~ 3.6 MeV ( 56 Co decay 
with an e-folding time 113.7 days). These line 7-rays 
are degraded along their paths through the SN ejecta by 
pair production, compton scattering, and photoelectric 
absorption (e.g., Casse & Lehoucq 1994 for a review). 
The absorbed energy is the source of the optical output 
of a SN in the case of Type Ib/c SNe (SNe Ib/c) (e.g., 
Maeda et al. 2003 for a review). 

By looking at properties of these 7- and X-rays, it is 
expected that we can directly extract information on the 
distribution and the amount of 56 Ni, which are closely 
related to the explosion mechanism. Such observations 
have been performed only for very nearby SN 1987A 
(Dotani et al. 1987; Sunyaev et al. 1987; Matz et al. 
1988). Attempts to detect the high emission from distant 
SNe (e.g., at ~ 10 Mpc) have not been successful because 
of limited observational instruments in the past (Matz & 
Share 1990; Lichti et al. 1994; Morris et al. 1997; Leis- 
ing et al. 1999). However, now that INTEGRAL and 
SUZAKU have been operated, and some new observato- 
ries are being planned (e.g., Takahashi et al. 2001), it 
is time to revisit theoretical prediction (e.g., McCray et 
al. 1987; Woosley et al. 1987; Shibazaki & Ebisuzaki 
1988; Kumagai et al. 1989) taking into account recent 
developments in SN research, i.e., multidimcnsionality 
(see e.g., Maeda & Nomoto 2003, Maeda et al. 2006ab, 
and references therein). 

Recently, multidimensional models for a supernova ex- 
plosion become more and more popular. Especially, bipo- 
lar explosion models have been intensively studied in 
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supernovae 

their explosion characteristics and nucleosynthetic fea- 
tures (e.g., Nagataki 2000; Maeda et al. 2002; Maeda & 
Nomoto 2003). There have been some observational evi- 
dence of such bipolar explosions - the HST image of SN 
II 1987A (Wang et al. 2002) and polarization measure- 
ments in SNe II/Ib/Ic (Wang et al. 2001; Kawabata et 
al. 2002; Leonard et al. 2002; Wang et al. 2003; Leonard 
et al. 2006). 

To judge a validity of a model, we have to compare 
model predictions with observations in details. For this 
purpose, we have been developing SupernovA Multidi- 
mensional RAdlation transfer code SAMURAI (Maeda 
et al. 2006ab; Maeda 2006; Tanaka et al. 2006). The 
code has been applied to very energetic supernova (hy- 
pernova) SN Ic 1998bw associated with a Gamma-Ray 
Burst. By comparing the model predictions with optical 
spectra and light curve of SN 1998bw in details, Maeda 
et al. (2006ab) concluded that SN 1998bw is the bipolar 
explosion viewed from the jet direction. 

In this paper, we present the model prediction for ra- 
dioactive decay 7- and X-rays from supernovae driven 
by bipolar jets. Details are presented in Maeda (2006) 
(see also Hungerford et al. 2003, 2005). Although the 
model has been applied to hypernovae, we believe that 
most of the model predictions apply even to usual SNe 
Ib/c which also show some evidence of asphericity (e.g., 
Wang et al. 2001). 

2. METHOD AND MODELS 

We use a fully time- and energy-dependent radiation 
transfer scheme to follow the high energy photons result- 
ing from 56 Ni/Co/Fe decay chain. This is done by using 
a Monte-Carlo method. Line 7-rays are produced ac- 
cording to the decay probabilities, then the photons are 
followed along their paths in which they suffer from pair 
creation, compton scattering, and photoelectric absorp- 
tion in expanding 3D SN ejecta. The detail is described 
in Maeda (2006). 

We adopt the bipolar model A from Maeda et al. 
(2006b), which successfully reproduces the optical obser- 
vations of hypernova SN 1998bw. The hydrodynamic and 
nucleosynthetic features have been examined in detail in 
Maeda et al. (2002). For comparison, we also exam- 
ine the spherical model S. The ejecta structure (density 
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Fig. 1.— Model A (Maeda et al. 2002) with E 51 = E/10 51 erg = 20. Left: Density distribution (g cm -3 ) on a logarithmic scale at 
t = 10 days. The ejecta are already freely expanding (so called homologous expansion or Hubble flow), so the distribution is shown in 
velocity space. Right: Mass fraction of 56 Ni on a liner scale. 



and 56 Ni) is shown in Figure 1. In Model A, 56 Ni is 
distributed preferentially along the jet axis because the 
materials in this direction experience high temperature 
(i.e., > 5 x 10 9 K) at the explosion. On the other hand, 
the density is enhanced in the equatorial plane (showing a 
disk- like structure), since only a weak shock wave pushes 
materials outward in the equatorial direction (Maeda et 
al. 2002; Maeda & Nomoto 2003). 

3. RESULTS 
3.1. Spectra 

Figure 2 shows the high energy spectra at t = 25 and 50 
days, which are obtained by our simulations. Hereafter, 
t is the time after the explosion. 

The aspherical model A has three unique features as 
compared to the spherical model S. 

(1) Large line-to-continuum ratio: Model A 
shows a large line-to-continuum ratio, especially 
early on. In Model A, the source of the high en- 
ergy emission, i.e., 56 Ni/Co, is abundantly present 
around the jet axis (z-direction) . In this direction, 
the density is small, thus the optical depth is small, 
as the consequence of the jet-induced explosion. 
Thus, few photons suffer from multiple compton 
scatterings which create the continuum in Model 
A. Since the situation does not change very much 
as time goes by, the line-to-continuum ratio does 
not show noticeable temporal evolution. 

Also, the line-to-continuum ratio is not sensitive 
to the viewing orientation. Both observers in the 
z- and r-direction look at the same region, i.e., the 
56 Ni/Co-rich jet. The line-to-continuum ratio is de- 
termined by local physical conditions in the emit- 
ting region, which is same irrespective of the ob- 
server's direction. Thus, it is nearly independent 
from the observer's direction. An exception is a 
line profile, which is discussed in §3.2. 

In Model S, a situation is different. The 56 Ni/Co- 
rich region is surrounded by massive envelopes in 



Model S. The line photons experience many scat- 
terings, creating the high level of the continuum. In 
Model S, the envelope becomes less and less optical 
thick due to the SN expansion, thus the continuum 
level becomes lower and lower as time goes by. 

(2) Large cut-off energy: The hard X-ray cut- 
off, which is created by photoelectric absorption, 
is at higher energy in Model A than Model S. In 
Model A, the last scattering point is close to the 
56 Ni-enhanced jet, thus the photoelectric absorp- 
tion takes place in the Fe-rich materials. Since the 
situation does not change as time goes by, the cut- 
off energy does not show temporal evolution. This 
behavior is again insensitive to the observer's di- 
rection. 

On the other hand, Model S shows temporal evo- 
lution. Early on, the last scattering takes place 
near the surface where O is the dominant element. 
Later on, the last scattering point moves deeper to 
the inner region where Fe is the dominant clement. 
Thus, the cut-off energy moves from lower to higher 
energy. 

(3) Early emergence: In Model A, escape of the 
high energy photons occurs earlier than in Model S 
if the kinetic energy of the expanding SN ejecta is 
the same. This is expected since 56 Ni is produced 
along the axis where the optical depth is low in 
Model A. 

3.2. Line Profiles 

Line profiles are also different for models A and S. As 
an example, Figure 3 shows the line profiles of 56 Co 1,238 
keV line. Other lines have similar properties, but the line 
blending should be taken into account (e.g., 56 Ni 812 keV 
and 56 Co 847 keV). 

Model A yields a blueshifted line if viewed from the jet 
(z-) direction. In this case, an observer views the opti- 
cally thin 56 Ni/Co-rich region from the jet (z-) direction, 
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Fig. 2. — Synthetic hard X- and 7-ray spectra (a) at t = 25 days and at (b) t = 50 days. 
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Fig. 3.— Profiles of the synthetic 56 Co 1,238 keV line (a) at t = 25 days and (b) at t = 50 days. 



and the emission from the other hemisphere is blocked 
by the optically thick central region. The central region 
continues to be optically thick at t < 100 — 200 days, 
thus the line keeps to be blueshifted. 

If Model A is viewed from the r-direction, then two 
56 Ni/Co-rich jets in both hemispheres can be seen. In 
this case, the emitting regions (the jets) are moving to the 
directions perpendicular to the observer. Thus, the line 
is centered at the rest wavelength when viewed from the 
r-direction. The situation does not change as time goes 
by, thus the line center does not evolve as a function of 
time in Model A (irrespective of the observer's direction). 

However, the situation is different for Model S. Early 
on, only the outer region of the 56 Ni/Co-rich region is 
seen, yielding the blueshifted line. As time goes by, the 
observer can look into the deeper region (because of the 
ejecta expansion and the density decrease), thus the line 



center continuously moves from the blue to the rest wave- 
length. 

4. CONCLUSIONS AND DISCUSSIONS 

We have shown in this paper that expected observed 
features in the high energy emission are highly depen- 
dent on the geometry of the explosion. Models A and S 
have quite different characteristics, thus these model can 
be distinguished by observations of the high energy emis- 
sion. Among the features, the line profile is very sensitive 
to the viewing orientation in Model A. Thus, if we can 
perform deep observations to detect the line profile, then 
even the viewing orientation can be distinguished. In ad- 
dition, temporal evolution is very different for models A 
and S. The expected features are summarized in Table 1. 

For the sensitivity of INTEGRAL (3<r in 10 6 sec), we 
estimate that the maximum distances within which the 
7-rays are detectable are ~ 500 kpc, 4 Mpc, and 7 Mpc 
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TABLE 1 
Characteristics 



Features Model S Model A - Z Model A - R 



Line-to-Cont a Small Large Large 

Cut-off E b Low High High 

Evolution of spectra c Large Small Small 

Line profile Evolved Blue Rest 

Evolution of Lines Blue to Rest Small Small 



a The line-to-continuum ratio, determined by compton scatterings. 
D The cut-off energy, determined by photoelectric absorptions. 
c The evolution of the Linc-to-continuum ratio and the cut-off energy. 
°~The line profile of isolated lines (especially at 1,238 kcV). 
e Temporal evolution of the line profile. 

for SNe II, canonical SNe Ib/c, and hypernovae, respec- 
tively. Since SN II 1987A, the nearest SNe for each type 
arc SN II 2004dj (~ 3 Mps), SN Ic 19941 (~ 7 Mpc), 



and the hypernova SN Ic 2002ap (~ 10 Mpc). Thus, the 
sensitivity of INTEGRAL is unfortunately not enough 
to detect these core-collapse SNe in the pace of 1 per 2 
decades, unless very nearby SNe like SN 1987A fortu- 
nately appears in the sky. 

If new generation hard X- and 7-ray telescopes are de- 
signed to achieve a sensitivity ~ 1 magnitude or two bet- 
ter than the current instrument (Takahashi ct al. 2001), 
then it is expected that the high energy emission from 
core-collapse SNe can be detected ~ 1 SN per year. 
This sensitivity will lead to a comprehensive study of the 
SN explosive physics through observing the high energy 
emission. 



The author is the JSPS Research Fellow. 
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